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Abstract 

Using a real-time formalism of equilibrium and nonequilibrium quantum-field theory, we de- 



^ ■ rive the reaction-rate formula for neutrino-conversion (z/ u') process and uu' annihilation 
process, which take place in a hot and dense QED plasma with background (anti)neutrinos 
out of equilibrium. Also derived is the formula for the inverse processes to the above ones. 
Using the hard-thermal-loop resummation scheme, we include the contribution from the 
coherent processes. The decay/production of a neutrino causes an evolution of its spatial 
distribution. A scheme for dealing with this evolution is presented. For the case of isotropic 
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neutrino distribution, numerical computation is carried out for the parameter region of 
type-II super-nova explosion. Defferential reaction rate exhibits characteristic peak struc- 
ture, which comes from the coherent processes. The contribution from the above processes 
to the decay or damping rate of a parent neutrino v is also studied. 

ll.lO.Wx, 13.10. +q, 13.15. +g, 13.35.Hb 



I. INTRODUCTION 

For the past two decades, properties of neutrinos in background media have attracted 
much interest (see, e.g., [|I|). Interactions of neutrinos with a thermal background cause a 
change in the properties of neutrinos. The dispersion relation is the quantity that describes 
this change. Notzold and Raffelt were the first who comprehensively analyzed the dis- 
persion relation of a neutrino in a thermal background, where, among others, the damping 
rate of an electron neutrino is computed, the rate which is related to the mean-free path and 
to the refractive index. Computation is performed by neglecting Pauli blocking effects and 
using the bare dispersion relation for participating electrons. Radiative decay of a massive 
neutrino has been analyzed in 0. 

It is by now well known |^,^ that, in hot and/or dense QED, the thermal propagators of 
a soft photon^ and a soft electron (positron) are drastically changed from those of respective 
bare counterparts. The salient feature is the appearance of the imaginary part in space- 
like-momentum region, which comes from Landau damping mechanism. The dispersion 
relations for soft photon and electron are also largely changed. An effective or improved 

soft particle is the particle that carries soft momentum (IQ'^I = 0(e\/T^~-i-/?))- Here — e 
is the electron charge and T (/u) denotes the temperature (chemical potential). 



2 



perturbation theory, called hard-thermal-loop (HTL) resummation scheme in which 

the above-mentioned effects are taken into account, is established just after the work [Q. 

Recently, in relation to the possible neutrino oscillation, neutrino-conversion processes 
have attracted much interest. For instance, in [0, the rate of neutrino conversion vl 
in a hot and dense QED (or electron-positron-photon) plasma is computed using the above- 
mentioned effective perturbation theory. The responsible interaction is the magnetic dipole 
neutrino-photon interaction, vl is assumed to be thermally distributed in the QED plasma. 

In this paper, we deal with a neutrino-conversion (z/ u') process and a uu' annihilation 
process, which take place in a hot and dense QED plasma. We also deal with inverse 
processes to them. The background (anti)neutrinos are out of equilibrium. In Sec. II, 
we derive the reaction-rate formula for these processes, on the basis of the nonequilibrium 
quantum-field theory (outlined in Appendix A) supplemented with the effective perturbation 
theory of hot and dense QED. Also presented is a procedure of determining spacetime 
evolution of (anti)neutrino distribution function. In Sec. Ill, for the purpose of illustration, 
numerical computation is carried out for the differential reaction rate for the case of isotropic 
neutrino distribution. We are interested in the temperature and baryon-number density 
regions of type-II supernova explosion (cf. rrie « T, fi « rriion- As in 

we neglect the effect of the ions. The contribution from the coherent processes exhibits 
a characteristic peak structure in energy distribution of a "decay neutrino" — v' for the 
neutrino-conversion process and v' for the vv' annihilation process. We then study the 
contribution to the damping rate of v. Numerical computation is carried out for the case 
where no neutrino exists in background. Section IV is devoted to discussions. In Appendix 
A, we present a formalism for dealing with neutrinos in out of equilibrium system, which is 
used in the main text. In Appendix B, we give the form for the self-energy-part-resummed 
photon propagator. Finally, in Appendix C, we compute the necessary quantities for the 
present study. 
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II. REACTION-RATE FORMULA AND DAMPING RATE OF A NEUTRINO 



A. Neutrino-conversion process 

We deal with the system that consists of a hot and dense QED plasma and neutrinos. We 
assume that the QED plasma is in thermal and chemical equilibrium, while the background 
neutrinos are not. Then, the whole system is out of equilibrium. A neutrino-conversion 
process of our concern is 

iy{K) + QED plasma i'\K') + anything. (1) 

We assume that v is massive, left-handed neutrino (with mass m), v' is massless and left- 
handed. The four-momenta K and K' are K = (E, k) with E = y/k"^ + w? and K' = (k', k'), 
respectively. The QED plasma is assumed to be at rest. The total reaction rate for the 
process (P contributes to the damping or "decay" rate T^i of a parent neutrino u. also 
receives a contribution from the relative process to (|l|), 

u{K) + v'{-K') + QED plasma anything, (2) 

where v' is an antiparticle of v' and —K' = {k', — k'). 
The region for various parameters of our interest is 

T, fi,E, k', rrie « Mw and m < 2me. (3) 

Here M]y (me) is the mass of W boson (electron). We expect that, in the region (^, the 
result is insensitive to m (< rrie)- As a matter of fact, we have confirmed that, at least in 
the range < m < 100 eV/c^ no visible change is detected in the results to be given in 
Sec. III. Thus, we set m = throughout in the sequel, so that K = {k, k). In the region (0), 
we may use the effective Lagrangian, which, after Fierz transformation, reads 

£eff = -2V2G Mi^Lv) {eif^Le) + h.c] . 

Here L = (1 — 75)/2, G is Fermi's constant, and underlined fields stand for the fields in the 
"weak-interaction basis". 
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For the QED sector, we employ the real-time formalism of equilibrium QED and, for 
the neutrino sector, we use the nonequilibrium real-time formalism of quantum-field theory 



Tni. is written as 10,11 



F^ = -^Tr 



2k 



dqo 



dqo ' 



(4) 
(5) 

(6) 



where U is the lepton mixing matrix and Sji^ is the (21)-element of the 2x2 matrix self- 
energy part T,^{X;K) of nonequilibrium u. (Here and in the following, denotes 2x2 
matrix in a "type-space" [P|,P|-pn|.) X stands for the spacetime coordinates of the center 
of the region, where the reaction takes place, the definition of which is given shortly. The 
0(G2) contribution to f^, f^^\ reads 



2k J (27r)4 
X Tr 



21 



(Q) 



LI^^^S^,';\X;K-Q)'ypL 



(7) 



where H^"^""^ is the (21) element of the one-loop matrix function 11^ whose (ij) element is 



n 



{W)al3 



(Q) = 8tG' 



Y+j 



i2ny 



Tr 



(^,J = 1,2). 



(8) 



The diagrammatic representation of f^;^'' is given in Fig. 1. In Eq. (|D, no summation is 
taken over i and j. In Eq. (|^), S'^i\X\P) is the (21) element of the matrix propagator 
Syi{X]P) of out of equilibrium i/', which is a Wigner transform of its configuration-space 
counterpart, Syi{x,y): 

d^P 



where X = {x+y)/2. S\f in Eq. (^ is the [ij] element of an equilibrium electron propagator 
matrix Sp. 



Neutrino sector 



From now on, we restrict our concern to the systems, in which the distributions of 
(anti)neutrinos are of quasi-uniform near equihbrium or of quasistationary. For such systems. 



the gradient approximation in the derivative expansion is sensible [jTO[ . A brief derivation of 
S,y'{X; P) to the gradient approximation is given in Appendix A: S^,{X; P) = Sly\X; P) + 
SD\X; P) with 5*^,?^ the leading term and S^)^ the nonleading or gradient term. Sly\X; P) 
reads 



5i?)(X;P) = L/4?(X;P), 

Aij(P) 
A;^(P)-A^(P) -A^(P) 



SS\X;P) 



-2mt{j>o)U{X-P)5{P^)A 



+ 1 



where 



fu'{X-P) 



P2±ze(po)0+' 
9ipo)N,,{X;po,p) 



(9) 



(10) 



A, 



V 



(11) 

(12) 



(l ±1 
±1 1 

Here Ni^/{X;pQ,p) [Np'{X; \po\, — p)] is the number-density function of u' [z/'] at the space- 
time point X with "energy" |po| and momentum p [— p]. The nonleading term Sly\X;P) 
is given in Eq. (|A2| ) with u' for u. In our formalism outlined in Appendix A, S^)^ is ©(G^) 
smaller than S^?-*. This is because S^^^ includes dx^fu'{X] P), which, as seen from Eq. ( |51D 
(with v' for I/), is proportional to since so is F^i. Then, for Sq^i ^ in Eq. (|^), we substitute 
the (21) element of the leading term S'f,?'' (Eq. (||)). 



Electron sector 

The electron propagator matrix Se{P), whose elements are in Eq. (|^), is of the form 
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5,(P)=M(P). (13) 

Here <Se{P) is given by the right-hand side (RHS) of Eq. (|^) provided that fyi{X;P) is 
replaced by an equihbrium distribution function 

/e(po) = ^(Po)iVe(Po) + ^(-P0)[1 - iV,+ (po)], 
1 



e(po-M)/T + 1 ' 

^--(^°) = e(bo|J/T + i - (14) 

We decompose njj^'*"'^((5), Eq. (H), into transverse (T), longitudinal (L), and vector- 
axial-vector interference {VA) parts. T- and L- parts are proportional to the thermal self- 
energy part of a photon in hot and dense QCD. Thus, we may write 



.G 



2 



-ze-^^OQ.nSf )(g)] , (15) 
nr(Q) = 2.e^^/-^(go-2po) 



g2 y (27 

x4)(p-Q)4;)(p), 



where 6^'^/'°" is a fully anti-symmetric pseudo tensor with e°^^^ = 1, S^j^ = 1,2) is the 
(ij) element of Se in Eq. (0), and "PxiL) standard projection operator onto transverse 

(longitudinal) mode, 

(q) = - E ^7"^/^(5^^' - g^g^) (g^ = gVg), 

Pf(Q)=^7"''-^-7'f(q). (16) 



^ij^^\Q) -^l- dH) (^^) element of the T/L-component of the thermal self-energy 

part of a photon. 

^21 (Q) {S = T, L, VA) are related to 0,^ the so-called Feynman self-energy part 

nSf)(g) ^ n(?(Q) + o{qo)ii[f{Q) + e{-qo)iii\Q) 

{S = T,L,VA), (17) 



through 

tU[%,^{Q) = 2[e{±qo)+nB{\qo\)]lran^F\Q) 

{S = T,L,VA). (18) 

Here nB^x) = l/(e^^'^ — 1) and "Im" means to take the imaginary part with Feynman 
prescription. In Appendix C, H^'* {S = T, L, VA) is computed within the approximation 
rrie = 0, which is a good approximation for a plasma with high temperature and/or density, 
rrie << T, /i. In Sec. Ill, we discuss to what extent the approximation mg = is good one. 
It is worth mentioning that Ilp'^\Q) vanishes for vanishing chemical potential, /i = 0. 

0{CP') decay-rate formula 

As seen from Eq. (0), S^I\X] K - Q) in Eq. (0) involves /^/(X; K-Q). Let us write 

fyi{X]K -Q) = /^/(X;/i;-go,k-q||,-q_L), 
q|| = ^(q-k)k, q± = q-q||. 

Then we define f,y' through 



qk 



J dqodq ^^^ q^f^^ (X; k-qo,k- q|[) g{qo, qy) 



(19) 



with Q any function of go and qy. 

Using the formulae displayed above, we obtain, after straightforward manipulation, for 
dVfldqQ to 0(G2) (ggg Eqs^ p ^nd (0)), 

df^}^ 1 1 



dqo TT^ k"^ 



[^(go) + ^B(|gol)]e(fc-go) 



X / dqq 
Jn 



qo , Q' 



1-/. X;.-,o, 1-f + ^ k 



G(i)(Q), (20) 



where 



kqo 



2 I -"iO TT 
Q \ 2 



(21) 



HvA{Q) = {qo-2k)Q\ 



The integration region TZ in Eq. (^OD is defined as (see Fig. 2) 

7^ = 7^l u 7^2, 

7^1 : |go| < g < 2A; - go, 

7^2 : \2k -qo\ <q< go- 
In the region TZi [7^2], = ^ ~ go ^ [^0 < 0]' ^^"^ then TZi [7I2] is the kinematically 
allowed region of the reaction (|l]) [(3)]. At first sight, at go = 0, Eq. ([20D seems to diverge 
at g = 0. Inspection of the formulae in Appendix C tells us, however, that this is not the 
case. 



B. Neutrino-production process 

We now turn to the inverse processes to and to (|]): 

u'iK') + QED plasma v{K) + anything, (22) 
QED plasma ^ v{K) + v\-K') + anything, (23) 

where K' = {k',k'), K = (fc,k), and —K' = (/c',— k). The process ( p2D is a production 
process of u due to the reaction of u' with constituents of the QED plasma and the process 
(p3|) is a uu' production process. The reaction-rate formula for these processes is written as 

Tp=\U:,Ue,,\'Tp, (24) 
9 



^ 2k 



df 

oo dqo 



In a similar manner as above, we obtain for the 0{G'^) contribution, 
dTW 1 1 

X J^dqqU, (X; k-qo,(l-^ + ^] k\ G«(g), 



(25) 
(26) 



(27) 



where fu' is as in Eq. (^) and G^^\Q) is as in Eq. (21). The diagram for f^^-' is the same 
as Fig. 1, provided that the two types of vertices are interchanged, 1^2. 



C. Contributions from a set of coherent processes 

According to the HTL-resummation scheme [^0, the integration region in Eqs. (|^), (^^, 
]), and (|27|) should be divided into hard-Q region {\Q^\ = 0{y/T^~+~jj^)) and the soft-Q 



region (|g^| = 0{ey/W+Ji^)). 

Hard-Q region: For Imllp (Q) (S = T,L,VA), expressions given in Appendix C are used. 
Soft-Q region: Observing the formulae in Appendix C, we see that, for e << 1, 

i/T/L(Q)Imn™(Q) ~ HT/L{Q)lmFT/L{Q) 

» HvAiQ)T^Z^\Q), (28) 



where ImFx/LiQ) is as in Eq. (|Cl| ) with Eqs. ( |C14|) and ( |C15|) in Appendix C. In the soft-Q 



region, there is an additional contribution: An inverse HTL-resummed photon propagator 
(cf. Eq. O) (*Af /^^(Q))"^ (= Q2 _ U^^/^\q) ~ Q2 _ Ft/l{Q)) is of the same order of 
magnitude as Ilp'^^\Q) (^ Ft/l{Q))- Thus, the diagram for as shown in Fig. 3 yields 
an equally important contribution. 

The characteristic scale of the hard region is im yr^T/?, and that of the soft region 



is ca/T^ + jj?. As a matter of fact, since e ~ 0.30, the hard region and the soft region are not 
sharply separated. Taking this fact into account, we compute the contribution from Fig. 3, 
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r(^\ without using the HTL-approximation (e << 1). The contribution is given by Eq. 
with the replacement MM, 



{W)af3 
21 



x(*A.,(g))^,n(r^(Q), 



(29) 



where II's are as in Eq. (H) and *A is as in Appendix B. Straightforward computation yields 



dr 



(2) 



dqo 



1 1 

[^(^o) + ^Bi\qo\)] e{k - go) 



X / dqq 
In 



(30) 



where 



Im 



P=T,L 



Q'-uPiQ) 



+HT{Q)q' 



(nr(Q) 



2Hva{Q) 



{VA). 



(31) 



g2-nP(g) g2_nP(g) 

The replacement ( pOf ) and Fig. 3 tell us |jll| that Eq. (|30D describes the differential rate for 
a set of processes, in which real and/or virtual photon(s) participate. It is to be noted that 
the (real) photons in the QED plasma are in thermal equilibrium. Then, the photon(s) does 
not come out of the plasma, so that, when the decay neutrino goes out from the plasma, it 
does not accompany photon(s). [In this relation, see [Q]. 

In a similar manner, we obtain, for the contribution from the processes (|2^) and (p3D, 

C 1 1 



dqo 



2 o2 J^2 



X 



[0{-qo) + nB{\qo\)]^iqo - k) 



(32) 



The diagram for this is the same as Fig. 3, provided that the type-1 vertex and the type-2 
vertex in Fig. 3 are interchanged. 

In the next section, we shall use the formulae displayed above for the whole g^-region. 
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D. Net decay rate 



The net decay rate r^*^* and the net production rate r^*^* are 

TT{X;k) = K{X;k)T,{X;k) 

-[1 - N,{X;k)]Tp{X;k) 
= -rf(X;k), (33) 

where Nj^lX^k.) is an (on mass-shell) distribution function of u. When i/, z/', and z/' are in 
thermal and chemical equilibrium, N,^, Nyi, and Nyi take similar form to Eq. (0). Using 
Eq. (|T^), one can show, in this case, that the detailed balance holds, V^^{X] k) V^^{k) = 
0. 



E. Procedure of determining (anti) neutrino distribution functions 

Here we study the neutrino or weak- interaction sector. As discussed above after Eq. (|l^), 
the nonleading piece of S„', Sl)\ is 0{G'^) smaller than the leading piece S^^^ and may be 
neglected. Higher-order corrections to Sgi^ Eq. (^, and to T^^i2 , Eq. (|25D , due to weak inter- 
action, are also 0{G^) smaller than the respective leading contributions analyzed above, so 
that they can be ignored. There is, however, an exception to this argument. Our formalism 
outlined in Appendix A accompanies, for each neutrino u, the evolution equation for f^, 
which is included in the neutrino propagator S,^{X; P): 

0xMX;P)L = F,{X;P), (34) 
F^iX; P) = z [(1 - /.)SS^ + . (35) 

Here S-J^ is the (ij) element of the self-energy-part matrix S,^ and is of 0{G'^). For com- 
puting Fi,{X;P), all the relevant contributions should be included, among which are the 
contributions analyzed above. (See Appendix A for more details.) 

We show in Appendix A that, on the mass-shell po = ±p, Eq. ( P^ ) turns out to the 
Boltzmann equation (|A7|) and its relatives, Eqs. (|A8|) and ([A9|). As a matter of fact, as far 
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as ©(G^) contributions are concerned, only /j, on the mass-shell is relevant because of the 
presence of 6{P^) in Eq. (^. Eq. ( plD describes the spacetime evolution of fu{X', P) under a 
given initial data fu{X^n^ X; P) at an initial time X° = X^^. Although the RHS of Eq. 
is of 0((j^), its effect cannot be ignored in general. This is because, in the process of solving 
Eq. for fu{X; P), integration over large spacetime scale is involved, so that the O(G^) 
effect is enhanced. 

In conclusion, neutrino self-energy part can be ignored everywhere but in in Eq. (^). 
In computing the reaction rates, as dealt with here, one should substitute the (anti)neutrino 
distribution functions, which are determined through Eq. (0) in a self-consistent manner. 
This applies to every neutrino. 



III. NUMERICAL COMPUTATION 

In this section we present a result of numerical computation of the "portion" of drd(p) / dqo 
and Td{p) that are independent of the z/' distribution function. We are interested in the type- 
II supernova environment, which is a QED plasma whose core temperature is T ~ 30 — 60 
MeV and electron chemical potential is /i ~ 350 MeV [Q. 



A. Differential reaction rates 



In general, the number-distribution function of z/' (z/'), N,^>(^pi-j{X; k'Q,k') {K' = K — Q), 
is anisotropic, and one should compute Eqs. (]20|), (pT]), (|30|), and ( p2D substituting N^i and 
Np/, which are to be determined self consistently. In this section, we restrict ourselves to 
the case of isotropic distribution, Ai',^/(p/)(/cQ, k') = A^jy'(p')(/i:o)- In this case, Eqs. (pO]) , (p^), 



and (|3^) may be written in the form. 



xg^'\qo,k,T) (^ = 1,2), 
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(2 = 1,2). 

Here 

g^^{q,XT)^e{k-q,)[e{q,)+nBm)] 



X 



/ dqqG^\Q) (^ = 1,2) 



with G^^\Q) and G^'^\Q) as in Eqs. (||) and (|3), respectively. We compute g^'\qo,k,T) 
{i = 1, 2) for various values for the parameters T and A;. For the chemical potential, unless 
otherwise stated, we take /i = 350 MeV 0]. 

In Figs. 4 - 8, we display the results of numerical computation for different values for k 
and T. The solid lines represent the total contributions Q = Q^^'' while the dot-dashed 

lines represent Q^^\ The figures "(a)" display Q and Q^^^ in the region qo < k [the region of 
the processes (0) and (P^)] and the figures "(b)" display Q and Q^^^ in the region qo > k [the 
region of the processes (0) and (|23|)]. Some observations are in order. 

• Figures 4-7 show the results for different values of k with T = 50 MeV. We see, as is 
expected, that for smaller incident-neutrino energy k, G^'^^G is larger. In the region 
of figures (b) [the region of the processes (H) and (^)], both G^^^ and G^"^^ {= G — G''^^) 
are positive. In the region of figures (a) [the region of the processes (0) and ( ^) ], 
except for the small region go ~ (or k' ~ k) in the case of relatively small k/T, 
G^'^'^ is negative. Referring to the reaction-rate formula [0, one can see what kind of 
physical processes are involved in df^^ /dq^ (Eq. (|30D) and df^^dq^ (Eq. (0)). As a 



p 

matter of fact, each of them involves a set of infinite number of coherent processes. A 
few examples of them that are involved in dT)i /dqo are 

z/-|-e— >e + 7 + i^', u + e^e + e + e'^ + h''. 



Figures 4-8 tell us that, for most regions displayed in figures (a), an infinite number 
of "interference contributions" is summed up to be negative, so that dTy/dqo and 
dT^ydqQ are negative. 
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• Both in figures (a) and (b), Q^^\qQ, k,T) exhibits peak structure. For figures (a) 
[go < k] , the peak is at go — or ^' = ^ ^ % — ^ and is more prominent for smaller 
incident energy k. The structure of figures (b) [go > k] may be understood as follows. 
In the hard-thermal-loop approximation [cf. Eq. (|28D ], *A^p^^\Q) ^ l/[<5^ — Ft/l{Q)] 
(see Eqs. ( p3|) and (|Cl|)). Then, as is well known [Q or as can be shown from Eq. ( P3|) 
with Eqs. and *A^^JI^^^'^\Q) in Eq. Q turns out to be of the form 

*A?;f ^^)(g) = 2t[9iTqo)+nBi\qo\)]ZT/Liq) 

x6{qo-ujT/Liq)) (go > g)- 

(36) 

The dispersion curves, go = ujriq) and go = ujiiq), are schematically shown in Fig. 2. 
Use of the actual *A^^^{Q) = 1/[Q'^ — Il^^^{Q)] results in the change of 5(go — i^T/L(g)) 
in Eq. (^) to the functions with finite width that are (sharply) peaked at go ~ oox/Liq)- 
Inspection of Fig. 2 with these observation in mind allows us to understand the struc- 
ture of figures (b). 

For the purpose of seeing the effect of the chemical potential fi, we display in Fig. 9 the 
result for {k,T,fi) = (10,50,0) MeV. [For the QED plasma in the early universe, /i ~ 0.] 
We see that Q^'^^ << G^^\ so that the peak structure is less prominent when compared to 
the case of /i = 350 MeV, Fig. 5. Q in the region go > is much larger than Q in the region 
go < k. 

Above computation is carried out neglecting the electron mass nie- Inclusion of the 
electron mass rrie causes a change in Il^p\Q) {S = T,L,VA) in the region \Q'^\ < O^ml). 
For the purpose of getting a measure to what extent the approximation mg = is good one, 
we perform all numerical computations by simply cutting off the region \Q'^\ < m^. This 
cutoff turns out to reduce Q^"^^ (i = 1,2). Dashed fines in Figs. 4-9 show the result of 
computation. In most regions of Figs. 4 - 9, no substantial reduction arises. Especially, for 
the region go > k, no sizable reduction arises and we do not display the results in figures 

15 



(b). For the region go < prominent reduction occurs only at go — 0, at which Q^^^ peaks. 
Larger reduction occurs for smaller k. 

B. Decay rate 

For computing the contributions to the decay or damping rate (= f ^^-^ + f^^-*) (cf. 
Eq. (ID) and to the production rate Tp (= T'^'> + F^^^) (see Eq. (p6D), knowledge for the 
distribution functions, N,^/ and Np/, is necessary. Furthermore, for computing the net decay 
rate Tf^^X; k), Eq. (0), knowledge for the distribution function N^{X] k) is necessary. 

Here we compute the damping rate of an incident z/ on a hot and dense QED plasma 
with no background neutrinos, = Ni,/ = Np = Npi = 0. Then, the process (H) is absent. 
Displayed in Figs. 10 and 11 are the total contribution f (= f^J^ + ff^) and the partial 
contribution f for fi = 350 MeV and, in respective order, T = 50 MeV and 20 MeV. 
Figures 12 and 13 show the result for fi = and, in respective order, T = 50 MeV and 20 
MeV. 

Cutting off the contribution from the region < ml does not result in sizable reduc- 
tion. 

Characteristic features: 

• For the range of k, T, and /i displayed in Figs. 10 - 13, the contribution from the 
soft-Q region, f^^"*, is not very large. 

• Figures 12 and 13 tell us that, for /i = 0, F^ is almost linear in k. As a matter of fact, 
Frf's in Figs. 12 and 13 are well parametrized as 

rd{k,T) = cG^k^+^T^-'' 

with (c,a) = (0.60,0.04). 
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IV. DISCUSSIONS 



On the basis of the formahsm, outhned in Appendix A, for deahng with nonequihbrium 
quantum-field systems, we have derived the reaction-rate formulae for neutrino-conversion 
process and its relatives, the processes which occur in the medium that consists of a hot 
and dense QED plasma and background neutrinos. The formalism involves the Boltzmann 
equation and two related ones, which describe spacetime evolution of (anti) neutrino distri- 
bution function. Illustrative computations of the differential reaction rate drd{p) / dqo is made 
in Sec. IIIA and of the total decay rate in Sec. IIIB. 

For relativistic particles dealt with here, the mean-free path / is related to the decay 
rate F through / = 1/F P], which, in turn, is related to the imaginary part of the refractive 
index Im[r;,] = {2lk)~^ = T /2k. Computation in Sec. IIIB shows that, in the range of Figs. 10 
and 11, order of magnitude of F^^ is 10^^^ ~ 10^^^ MeV. Then, we see from Eq. (^) that 
/ ~ [0.2 ~ 200]/|f/*j,f/e!/'P m, which is much less than the core size of the type-II supernova. 
This means that, when applying to the actual supernova, vv' as well as vv' production 
processes are also important and, through these processes, (anti)neutrinos are produced. 
Thus, one has to perform an analysis by the full use of the formalism in Appendix A. 
The evolution of the neutrino- and antineutrino-distribution functions should be dealt with 
through the Boltzmann equation and its relatives. Once the initial distribution function of 
(anti)neutrinos are given, the framework presented here allows one to determine hereafter of 
the system. Concrete numerical analysis along this line is outside of the scope of the present 
paper. 
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APPENDIX A: NONEQUILIBRIUM NEUTRINO PROPAGATOR AND 
"HEALTHY" PERTURBATION THEORY 



A perturbative framework for dealing with massless Dirac fermions in a nonequilibrium 



quantum-field system is presented in |12|. Extension of the framework to the case of mass- 
less left-handed neutrino vl is straightforward. Here we briefly describe somewhat simpler 
framework for z/^, which is sufficient for the present purpose. We employ the derivative 
expansion and use the gradient approximation throughout. 



1. Free propagator 



We start with introducing a standard form |[13|JT^ for nonequilibrim matrix propagator 
of a massless left-handed neutrino vl. 



BR{x,y) 



X S^Aiu - v)BR{v,y), 
S^{x-y) ~fu{x,y) ^ 

^5^{x-y) 6^{x - y) - fy{x,y) ^ 



(Al) 



5^{x - y) - f^{x,y) -f^{x,y) 



\ 



6'^{x - y) 5'^{x -y) j 

Sra = diag {Sr, -Sa) ■ 



Here denotes a 2 x 2 matrix in a "type space." In a 4 x 4 Dirac-matrix space, Bl and 
Br are the unit matrices. The Fourier transform of Sr(^a) reads 

Sr(A){P) = Lf^RiA){P) {L = {1- 75)/2), 

where /S.r(^a){P) is as in Eq. ([T0|) and fu{x,y) is the inverse Wigner transform of 

UX;P) = 9ipo)N,{X;po,p) 

+9(-p,)[l-N,{X; bo|,-p)] 



with X = {x + y)/2. Here N,y (Np) is the number-density function of u (z/). Computation 
of Eq. ( |A1| ) to the gradient approximation yields 



-iP-{x—y) Q 



X + y 



(27r)4 ""V 2 

S,iX;P) = Sl'^\X;P) + Si'\X;P). 



The leading term S^^^ is given by Eq. (P) with u for u' and the nonleading term S^^^ reads 



p2 



Qp2 p2 



(A2) 



where is as in Eq. ([12D and P denotes principal part. 



2. Free and counter actions 

The propagator introduced above is an inverse of the kernel of a free action Ao, which 
we now find. Applying 0^ to Su{x,y) in Eq. (|A1|) , we find the form of Aq, to the gradient 
approximation: 



A, 



d*x^{x)i$m{y) ~ Ac, 



(A3) 
(A4) 



d'^x f d^y 



(27r)4y (27r)4 

X {z0, + i0y)Lf,{x,y)A_^{y), 



(A5) 

where A_ is as in Eq. (|12D. In Eq. ([A^), the subscripts of the field denote the type of field 
in real-time formalism [T^. From Eq. (|A^), we see that the action of the theory turns out 
to be 

A = Ao + A' + Ac, 

where A' includes other fields than ul and ul and the interactions between fields. It is to 
be noted that the counter action Ac appears. 



19 



3. Self-energy part and "healthy" perturbation theory 



The self-energy-part matrix of v reads 

%{X- P) = P) - i$xf{X; P)LA_, 

where S(,'°°p) comes from loop diagrams and the second term on the RHS comes from the 
counter action Af.. Computation of a S-inserted propagator to the leading order of derivative 
expansion yields 

S,{X-P)t,{X-P)S,{X-P) 
+ .... 

Here = f^{X;P) and contains the terms with Sj^Tj^^^Sr and with Sa^a^Sa, where 
— ^11 + ^12(21) is retarded (advanced) self-energy part [|1^. Observing SrSa oc 
1/ [(P^ + iO^)(P^ — iO^)], we see that SrSa possesses pinch singularities at po = ±p in a 
complex po-plane. Then, by demanding^ 

0xfu{X;P)L = F,{X;P), 



F,{x-p) = i (1 - /,)ssj + /,sri 



(A6) 



we attain a pinch-singularity free perturbation theory. 



4. Boltzmann equation and its relative equations 

After multiplying from the left of Eq. ([A6|) , we take a trace and go on the mass-shell 
Po = =tp- Referring to Eqs. @ and (^) together with their antiparticle counterparts, we 
obtain, with obvious notation. 



**As a matter of fact, demanding Eq. ( |A6| ) to hold in any region(s) of P will do, as far as po 
and Po = ~P are within that region(s). 
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= +p 



(5xo + v-Vx)iV±(X;p,p) 

= (l-iV±(X;p,p))rW(X;p) 

-iV±(X;p,p)rf)(X;p). (AT) 



Here v = p/p and "+" ["— "] stands for v [z/]. As can be seen from Eq. (|33D , the RHS is the 
net production rate oivjv and Eq. ( |A7|) is nothing but the Boltzmann equation. 

Similarly, multiplications of L^^ and and of L^^ (= L7 — L(7 ■ p)p) to Eq. (|A6|) yields, 
in respective order, 

axoAr±(X;p,p) = ±iTr[L7°FJ , (A8) 

7x ■ VxAr±(X; p) = ±i Tr [L7±F.],„=±,,p^±p • 

(A9) 

If (initial) distributions are spatially isotropic, the RHS of Eq. ( |A^ ) vanishes. 

APPENDIX B: SELF-ENERGY-PART RESUMMED PROPAGATOR OF A 

PHOTON 

Elements of the self-energy-part-resummed photon propagator matrix (in Landau gauge) 
reads [|,§,| 

*Ay(Q) = -Pf (q) *A?^)(Q) - n\Q) *A?'^(Q) 

(^,J = 1,2), (Bl) 

where V^^{c[) and Vl^{Q) are as in Eq. (|T6|) and 

*a(!;1(q) = -(*a?;1(q))* 

= *Ag^/^)(Q) + 2mB(|go|) Im*AS.^/^)(g), 
*Ag;i/(^^)(g) = 2z[^^(Tgo) + ns{\qo\)] Im*AiT/^)(Q), 

(B2) 

*Af /^)(Q) = . (B3) 
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Up ' is computed to the one- loop order in Appendix C. 



APPENDIX C: SELF-ENERGY PART H^f^ (5* = T,L,VA) 



Here we compute the lowest-order contribution to Ilp{Q) {S = T, L, VA)), Eq. ([T7|) . We 
are interested in the high-T and large-/i region T,fi » rrie, and then we ignore rrie- The 
effect of me (7^ 0) is discussed in Sec. III. 



Computation of n'^PiQ) {S = T, L, VA) 



We decompose Hp (Q) and Hp {Q) into three parts, 



/^)(Q) = F^7^(Q) + Ft/l{Q) + Gt/l{Q). 



(Cl) 



fI^Ij^ stands for the vacuum contribution and Ft/l stand for the contributions that domi- 
nate in the soft-Q region, the latter contributions which are called hard thermal loop 
Incidentally, n^^'*((5) has no hard thermal loop. 

Straightforward computation of Eq. (|l^) (cf. Eqs. (H), (|15D, and ([T6|)) using Eqs. (|13D 
and dH) yields 



FP{Q)=Fr{Q) = -gQ' 



FAQ) 

Fl{Q) 



GLiQ) 



a 



3 2 10 
2 ' q 



go _ qo + q 

q 2g2 go - g 



-3 m — 

1 gZ 



a 



Gt{Q) = —Q' 

TT 



TT q^ 



1 _ ^ In ^0 + ^ 
2q qo-q 



Q^h+Aih-qoh 



(C2) 
(C3) 
(C4) 

(C5) 
(C6) 



and 



n 



(VA) 



/igo In ggo^s + 2g/i 

go - g 



(C7) 
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In obtaining the vacuum contribution (|C2|) , we have used the MS scheme^^ and fir is the 
renormahzation scale, for which we choose fir = a/T^ + fi^- Incidentally, the vacuum part of 
n^'^'' vanishes. In the above equations, a = e^/An, = e^(T^ + 3;U^/7r^)/9 is the thermal 
mass of an electron, and 



with 



ap[n+[p) + n^[p)\ in 



2q Jo 



TC 



1 f°° ( 

~2qJo ^^^^ l'^+(P'> + ^-(P)] 1^ ( 



1 r°° ( 



h = dpp [nM - n-{p)] In (^^^ 



y^^dp[n4p)-n4p)] In 



n±{p) = l/(e(^^^)/^ + 1), 

Lp„ = qQ + pq + 2ap (p,cr = ±). 



(C8) 



(C9) 

(CIO) 

(Cll) 
(C12) 



It is straightforward to obtain 



Imni.^/^)(g) = - [nS.^/^)(go(l + ^e), q) - c.c. 



lmFT{Q)=e{-Q^)—nlQ 



2\ 371" 2 r^2 ml 



(C13) 

(CM) 



^We have adopted a convention that Dirac gamma matrices are 4x4 matrices in D-dimensional 
spacetime. 
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where 



Im72 = -— / rfpp^[n+(p) 

Zq Jqi 

Im/a = -—e{qo) / cipp[n+(p) + 



■-e(go)^(-g') 



Jo 

Im/i = -— / dpp [n+(p) - n_(p)] , 

ZQ Jnj 



2q 

Im/3 = ^e{qo)e{-Q^)fi 

-^e{qo)T[F^+-F_ 
2q 



_ \Qo\ + Q 



Qu = — ^ — , qi = 



2 ' 2 ' 



Note that U.^p'^\Q) vanishes for /j, — 0. 
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FIGURES 

FIG. 1. Diagram for f^^^- "1" ("2") at the vertex on the left-side (right-side) denotes the 
type of vertex in real-time equilibrium/nonequilibrium quantum field theory. The oval loop is an 
electron loop. 

FIG. 2. Integration region TZ = TZi U TZ2- The dashed line with T (L) shows the disper- 
sion relation for the transverse (longitudinal) mode in the hard-thermal-loop resummed photon 
propagator. 

~ (2) 

FIG. 3. Diagram for F^ . The hard-thermal- loop resummed effective photon propagator is 

indicated by a blob. "1", "2", "z", and "j" on the vertexes denote the type of vertex. The oval 

loops are electron loops. 

FIG. 4. Plots of g and ^(i) vs go at T = 50 MeV, fi = 350 MeV, and ^ = 2 MeV. Figure (a) 

corresponds to the processes (|l]) and (22), and Fig. (b) corresponds to the processes d) and dH). 

FIG. 5. Same as in Fig. 4 but for E = 10 MeV. 

FIG. 6. Same as in Fig. 4 but for ^ = 20 MeV. 

FIG. 7. Same as in Fig. 4 but for E = 50 MeV. 

FIG. 8. Same as in Fig. 4 but for T = 20 MeV and E = 10 MeV. 

FIG. 9. Same as in Fig. 4 but for ^ = and E = 10 MeV. 

FIG. 10. Plots of frf and f^^^ vs at T = 50 MeV and fi = 350 MeV. 

FIG. 11. Same as in Fig. 10 but for T = 20 MeV. 

FIG. 12. Plots of f and f^^^ vs at T = 50 MeV and fi = 0. 

FIG. 13. Same as in Fig. 12 but for T = 20. 
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